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Abstract 

A Shell-model-Like APproach (SLAP) suggested to treat the pairing correlations in relativistic mean 
field theory is introduced, in which the occupancies thus obtained having been iterated back into the 
densities. The formalism and numerical techniques are given in detail. As examples, the ground 
state properties and low-lying excited states for Ne isotopes are studied. The results thus obtained 
are compared with the data available. The binding energies, the odd-even staggering, as well as the 
tendency for the change of the shapes in Ne isotopes are correctly reproduced. 
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I. INTRODUCTION 



The relativistic mean field (RMF) theory is one of the most successful microscopic models 
111]. During the past two decades, it has received wide attention due to its success in describing 

nn 

many nuclear phenomena for the stable nuclei |2J, |3J as well as nuclei even far from stability 
It has been shown that the relativistic Brueckner theory can reproduce better the nuclear 
saturation properties (the Coester line) in nuclear matter |6j , present a new explanation for the 
identical bands in superdeformed nuclei []| and the neutron halo ^| , predict a new phenomenon 
- giant neutron halos in heavy nuclei close to the neutron drip line jjj], give naturally the 
spin-orbit potential, the origin of the pseudospin symmetry llfll [ll| as a relativistic symmetry 
Il^i ll^ and spin symmetry in the anti-nucleon spectrum |lf| , and present good description 
for the magnetic rotation and the collective multipole excitations Q], etc. 

Pairing correlations play an essential role in nuclear properties such as, binding energies, odd- 
even effects, single particle orbit occupation, electric-magnetic transition, low-lying collective 
modes, moment of inertia, halo phenomena, etc. In order to be realistic and to describe the open 
shell nuclei, the RMF theory must be supplemented with the proper treatment of the pairing 
correlations. This aspect becomes all the more important when considering the continuum 
effects, which are crucial for the description of drip line nuclei. So far the most commonly used 
methods, the Bardeen-Cooper-Schrieffer (BCS) approximation and Bogoliubov transformation, 



have become standard in the nuclear physics literatures 



18l . LL9f , even for exotic nuclei where 
the Bo g o liub ov —.nation is reeled to bejpn^ooriinate space flflHQorthe 

3- 



introduction of the resonant states is necessary [22], |23|, |24 , 

As the approximate product of the quasi-particle wave functions in quasi-particle formalism 
breaks the gauge symmetries connected with the particle numbers, the projection methods are 



introduced to restore the broken symmetries 



26 



There exists a vast amount of literature on 
and the references therein). The variation 



such projection methods (see, for instance, Ref. 
after projection (VAP) [28j] is the appropriate tool that fulfills the variation principle and 
provides a self-consistent description of fluctuations going beyond mean-field. Although, the 
method of variation after projection has been known long time ago, the numerical solution of 
the corresponding variation equations is relatively complicated. A fully self-consistent variation 
after exact projection has been carried out so far only in a limited number of cases (27, 29]. Exact 
projection within full HFB-theory is possible by a search of the minimum in the projected-energy 
surface by gradient methods [3^ . However these methods are numerically very complicated and 
have been applied only to the case of number-projection in restricted spaces and for separable 
forces The problem whether simple HFB-like equations can be obtained with the projected 



energy functional has been recently addressed arid it has been shown that the variation of an 
arbitrary energy functional, which is completely expressed in terms of the HFB density matrix 
p and the pairing tensor k, results in the HFB like equations with modified expressions for 
the pairing potential and the Hartree-Fock field j^. This would allow the usage of numerical 
algorithms for diagonalization of matrices as is done in the ordinary HFB theory. 

Meanwhile, the defects including the violation of the particle number conservation and im- 



probable treatment for the Pauli blocking effects |34j, |35|, |36|, in both the BCS approxi- 
mation and Bogoliubov transformation for finite fermion system like nuclei can be avoided in 
the Shell- model-Like APproach (SLAP) [3^], which was originally refereed as Particle Number 
Conserving (PNC) method there. 

Compared with the complicated particle number projection technique, the SLAP proposed 



in the 1980s 



34j ] avoids all the difficulties encountered in BCS and Bogoliubov approximation 



and takes into account the Pauli blocking effects strictly by diagonalizing the pairing Hamil- 
tonian directly in a reasonable configuration space. It has been found that the number of the 
configurations with important weights in the low-lying solutions is quite limited, only a few 
lowest-energy configurations contribute significantly 40]. This conclusion has been used as the 
guideline of a many-body configuration truncation by retaining a certain number of the lowest 
energy states. Extensive studies and discussion on the validity of the truncated many-body 
spaces as well as its application can be found in Refs. |44j and references therein. 

The application of the SLAP to solve the nuclear mean-field plus pairing Hamiltonian problem 
with a realistic deformed Woods-Saxon single-particle potential has been reported in Ref. j^. 
The optimal basis construction has been discussed and the stability of the final result with 
respect to the basis cutoff has been illustrated. In particular the presence of the low-lying se- 
niority s = solutions, that are usually poorly described in terms of the standard BCS or HFB 
approximation, has been found to play a role in the interpretation of the spectra of rotating 
nuclei. 

As SLAP is based on the direct diagonalization of the pairing Hamiltonian in the multi- 
particle configuration space, it has been shown to be more accurate than the BCS calculation 
as compared with the exact solution at least for the small systems j^. Furthermore, the 
blocking effects are taken into account automatically and both odd-A and even-even nuclei 
can be treated on the same footing. In addition, the excitation spectrum can be obtained 
conveniently in the present algorithm. 

In this paper, SLAP is suggested to treat the pairing correlations in relativistic mean field 
theory, in which the occupancies thus obtained having been iterated back into the densities. 
A brief formalism is presented in section II. The numerical details and its application for the 



ground state properties and low-lying excited spectrum for Ne isotopes as well as the comparison 
with data are given in section III. Finally, a brief summary is given in section IV. 

II. FORMALISM 

A. Relativistic Mean Field Theory 

In the framework of RMF theory, the nuclear effective interaction is usually described by the 
exchange of three mesons: the scalar meson a, which mediates the medium-range attraction 
between the nucleons, the vector meson c<j m , which mediates the short-range repulsion, and the 
isovector- vector meson p 5 ", which provides the isospin dependence of the nuclear force. The 
effective lagrangian density is as the following: 



where the field tensors of the vector mesons and the electromagnetic field take the forms: 




- T^O" 4 - g<r1p(7ll> 



(1) 



and other symbols have their usual meanings [l|, |2|, |3( . 

The classical variation principle leads to the Dirac equation, 





. V + V(r) +(3[M + S(r)}} fa = e t fa 



(3) 



for the nucleon spinors and the Klein-Gordon equations, 



[-A + ml]a(r) = -g a p s (r) - g 2 cr 2 (r) - g 3 a 3 (r), 
[-A + m>"(r) = gj^r) - g^uju^r), 
[-A + mj] p"(r) = gj^r), 



(4) 



-AA»(t) = e#(r), 



for the mesons, where 




(5) 



are respectively the vector and scalar potentials and the source terms for the mesons and the 
photons are 
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In Eqs.(6), the summations are taken over the valence nucleons only, i.e., no-sea approximation 
is adopted. The coupled equations (3) and (4) are non-linear quantum field equations, and 
their solutions are very complicated. Therefore, the mean field approximation is generally 
used to solve Eqs.(3) and (4): i.e., the meson field operators in Eq.(4) are replaced by their 
corresponding expectation values, and the nucleons are considered to move independently in the 
classical meson fields. In such a way the coupled equations (3) and (4) can be self-consistently 
solved by iteration. 

The symmetries of the system simplify the calculations considerably. In the system consid- 
ered in this work, as there exists time reversal symmetry, there are no currents in the nucleus 
and the spatial vector components of oo, pand A vanish. This leaves only the time-like compo- 
nents u!°, p° and A . Charge conservation guarantees that only the 3-component of the isovector 
p survives. For simplicity, we denote it in the following simply by p°. 

For axially deformed nuclei, i.e., the systems which have rotational symmetry around a 
symmetrical axis (assumed to be the z-axis), it is useful to work with cylindrical coordinates: 
x = r_i_ cos (f,y = r± sin ip and z. The potential of the nucleon and the sources of meson fields 
depend only on the coordinates r± and z. The Dirac spinor ipi for the nucleon with the index i 
is characterized by the quantum numbers O i5 7i~j and tj, where Qi = mi- + m s . is the eigenvalue 
of the angular momentum operator J z , 7Tj is the parity and ti is the isospin, i.e., the Dirac 
spinor ^ has the form [rs| . 

/ f+(z, r _ L ) e i (^- 1 /2)v> \ 

tg+iz^r^e 1 ^- 1 / 2 ^ 
\ig~(z,r ± )e i ^ +1/2 ^ J 

For each solution with positive ipi, there is the time-reversed solution with the same energy, 
ipl = Tipi, in which the time reversal operator T = —ia y K and K is the complex conjugation. 
For nuclei with time reversal symmetry, the contributions to the densities of the time reversal 




(7) 



states, i and i, are identical. Therefore, the densities can be represented as 



Ps,v 
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(8) 



and, in a similar way, p 3 and p c . The sum here runs over only states with positive Qj. The 
occupation number for state i is represented by rij. These densities serve as the sources for the 
fields cf) = a,uj°,p° and A , which are determined by the Klein-Gordon equation in cylindrical 
coordinates. 

To solve the equations (3) and (4), the spinors fi(z,r±) and gf(z,r±) can be expanded in 
terms of the eigenfunctions of a deformed axially symmetric oscillator potential |l8| or Woods- 
Saxon potential an d the solution of the problem is transformed into a diagonalization of a 
Hermitian matrix. The details can be found in Ref . [rsl. 4s| . 

The total energy of the system is: 
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where ^nucieon is the sum of the energy for nucleon e i} E a , E u , E p , and E c are the contributions 
of the meson fields and the Coulomb field, and -Ecm is the center of mass correction. 



B. SLAP for the Pairing Correlations 

Based on the single-particle levels and wavefunctions obtained from the RMF theory, 
the SLAP can be adopted to treat the nuclear pairing correlations, here after refereed as 
RMF+SLAP. The Hamiltonian reads: 

H = i/g.p. + impair, (10) 

where, H sp = ^2e v a^a u and H pSuir = —G a ^ a p a & a v, with e v the single particle energy 

v fi,v>0 

obtained from the RMF theory and i^ pa ir the pairing Hamiltonian with average strength G 



simplified from the quantization of the meson fields, v is the notation of the level with the 
quantum numbers (Q, ir, t) for axially deformed nuclei, and v is the time-reversal state of v. 

The multi-particle configurations (MPC) are constructed as bases to diagonalize the Hamil- 
tonian in Eq.(lO). For a even N = 2n particle system, the MPC are constructed as the 
following: 

(a) The fully paired configuration (seniority s = 0): 

\P1P1 ■ ■ ■ PnPn) = a+ a± ■ ■ ■ a + pn a+ n |0) , (11) 

(b) The configuration with two unpaired particles (seniority s — 2): 

\pvp\p\ ■ ■■pn-ipn-i) = a+a+a^a^ ■ ■ ■ o+_ 1 o±_ 1 |0) , (12) 



and so on |34| . The MPC for a odd N = 2n + 1 particle system can be constructed similarly 
For the axially deformed nuclei, as the parity it, seniority s, and eigenvalue K of the third 
component for the total angular momentum operator J z are good quantum numbers, the MPC 
space of a even particle system can be written as the following: 

MPC space = (s = 0, K = + ) 

© [s = 2, K = 0+) © (s = 2, K = 1 + ) © (s = 2, K = 2+) © • • • 

© ( a = 4, K = 0+) © (s = 4, K = 1+) © (s = 4, K = 2+) © • • • 

©•••. (13) 

In realistic calculation, the MPC space has to be truncated and an energy cutoff E c is in- 
troduced. Only the configurations with energies E{ — Eq < E c are used to diagonalize the 
Hamiltonian H in Eq.(10), where E^ and Eq are the energies of the ith configuration and the 
lowest configuration, respectively. The cutoff in the single particle states is implicitly included 
in the energy cutoff E c . 

The corresponding nuclear wavefunction can be expanded as 

PI,— >Pn P", v PU— >Pn-l 

where (3 — (ground state), 1, 2, 3, ■ • • (excited states). 
The occupation probability of the ith-level for state f3 is 

2 



E K-,p^ +E E K^L,,i +•■■' i=i,2,3,.... as) 



p\,---,p n -l pi,---,p n -2 



Replacing the occupation number m in Eq.(8) by nf~ , the sources in Eqs.(4) are obtained 
and the occupancies thus obtained having been iterated back into the densities, which give the 
new meson fields and new electromagnetic field. These new fields then are used to calculate 
the vector and scalar potentials in Eqs.(5). With these new potentials, the Dirac equation 
Eq.(3) is solved again. These processes should be repeated until the results converge to the 
given precision. As the pairing correlations are taken into account, the pairing energy -E pa ir = 
(ip 13 1 impair | ip 13 ) should be added to the total energy in Eq. Q. 



III. RESULTS AND DISCUSSIONS 



We have performed the RMF+SLAP calculations for the Ne isotopes with the usual effec- 
tive interactions, e.g., PKlQ], NLSH 

H, NL3 Q, and TMaQ- 

As the conclusions do not 

depend on the effective interactions used here, only the results with the effective interaction 
PK1 are presented in detail. Similar investigation has been done with the RMF+BCS approach 
in Ref. [^fj. For the axially deformed nuclei, the parity 7r and the eigenvalue K of the third 
component of total angular momentum operator J z are good quantum numbers. In the numer- 
ical calculations, the Dirac equation Eq.(3) for the nucleons and the Klein-Gordon equations 
Eqs.(4) for the mesons and the photon are solved by expansion in the harmonic oscillator basis 
with 14 oscillator shells for both the fermion fields and the boson fields. The oscillator fre- 
quency of the harmonic oscillator basis is fixed as ftwo = 41A -1 / 3 MeV and the deformation of 
harmonic oscillator basis (3q is reasonably set to obtain the lowest energy. 

Taking 24 Ne as an example, the detailed procedures of SLAP for the pairing correlations are: 

• From the solution of the Dirac equation, the neutron and proton single-particle levels 
are obtained, as schematically shown in Fig. 1, where the dashed line represents the last 
occupied level in the filling approximation. 

• With the single-particle levels thus obtained, the MPC are constructed according to the 
Eqs.(ll) and (12) ( see also Ref. Q] ) with given truncation energy E c . The truncation 
energy E c is chosen in such a way that no configurations with important weights are 
omitted in the SLAP calculations. In the following calculations on the Ne isotopes, the 
truncation energy is reasonably fixed to E c = 40 MeV. With this truncation energy, 
all the configurations with weight above 3.7 x 10~ 6 have been included for 24 Ne, while 



configurations with weight above 0.032% included 
the configuration truncation can be found in Refs. 
references therein. 



or Ne. More detailed discussion on 
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42J, 3, Q and the 



As the parity n, seniority s, and K are good quantum numbers, the Hamiltonian H 
in Eq.(lO) is block diagonal in the MPC space Eq.(13). In every subspace, the single 
particle Hamiltonian i? s . p . is diagonal. The non-zero off-diagonal matrix elements come 
from impair for configurations with a difference of one pair particles in the particle-pair 



occupation orbital. The pairing strength G T =C T /A MeV (r = p,n) is adop 



ted for all 



.[51]. 



Here 



the Ne isotopes, as normally done in RMF+BCS calculations, e.g., see Ref. 
for truncation energy E c = 40 MeV, C T = 10 is fixed by reproducing the experimental 
odd-even mass difference of 24-26 Ne. 

Diagonalizing the matrix of the Hamiltonian H in Eq.(10) obtained above, the solutions 
for the ground state and the low-lying excited states are obtained. The occupation proba- 
bilities nf of the single-particle level i in state (3 can be calculated from Eq.(15). Replacing 
the occupation number by n^~° in Eq.(8), the densities in Eqs.(6) can be obtained, 
which give the new meson fields and new electromagnetic field. These new fields can be 
used to calculate the vector and scalar potentials in Eqs.(5). With these new vector and 
scalar potentials, the Dirac equation Eq.(3) is solved again. 

The above processes should be repeated until the results converge to the given precision. 
The calculated results obtained in such a way will be refereed as RMF+SLAP in the 
following. 



nucleon 



As the pairing correlations are taken into account, for state /3, the energy for nucleons E-, 
in the total energy of system Eq.(9) should be modified as: 

rr -STri 13 * 

-'-'nucleoli / "'j fc i7 

i 

and the pairing energy -E pa ir = ( - /3 |-ff P air|V' /3 ) should be added to the total energy. 

The results for the ground state of 24 Ne from the RMF+SLAP* and RMF+SLAP calcu- 
lations with PK1 are shown in Table I, including the binding energy E, binding energy per 
nucleon E/A, the rms radii for neutron, proton, and matter, R n , R p , and R m , and the quadru- 
ple deformation for neutron, proton and matter, fi 2n , /3 2 p and /? 2m - The difference between 
these two calculations lies in that in RMF+SLAP*, the occupancies obtained by SLAP is not 
iterated back into the densities. While in RMF+SLAP, the occupancies obtained by SLAP is 
self-consistently iterated back into the densities. The corresponding calculations and the data 
available 0, are included for comparisons. From Table I, it can be seen that as pairing 
correlations are taken into account, both the RMF+SLAP* and RMF+SLAP calculations give 
more bound for the total binding energy and larger rms radii due to the contribution from 



those loosely bound levels, which is similar as that obtained from the Bogoliubov transforma- 

n 

tion, i.e., the pairing correlations can increase the nuclear rms radiusjqj. However, compared 
with RMF+SLAP*, the binding energy from RMF+SLAP calculations is less. Furthermore, 
the deformation of nucleus reduced by the pairing correlations is realized only in RMF+SLAP, 
which clearly demonstrated the necessity to iterate the occupancies obtained by SLAP self- 
consistently back into the densities. 

As the Pauli blocking effects can be strictly treated in SLAP, a self-consistent description 
for even-even, odd-A, and odd-odd nuclei can be obtained in the RMF+SLAP calculations. In 
Fig. 2, the occupation probabilities for the neutron levels are displayed for 24 Ne in the left panel 
and for 25 Ne in the right panel, where the occupation probabilities of levels above [202]5/2 + for 
24 Ne and those of levels above [200]l/2 + for 25 Ne are multiplied by a factor of 10 for clarity. 
For 24 Ne, if without pairing correlation, the levels above [202]5/2 + will be unoccupied and the 
others are fully occupied. In the RMF+SLAP calculations, the levels above [202]5/2 + have been 
occupied partly, while the occupation probabilities for levels below Fermi surface are decreased 
accordingly. For 25 Ne, the [200]l/2 + is the blocking level and is occupied by a single nucleon 
in the RMF+SLAP calculations. The occupations for the levels beyond [200]l/2 + have been 
reduced due to the blocking effects. For 24 Ne, both RMF+SLAP* and RMF+SLAP calculations 
have been done and it turns out that their difference in the occupation probabilities can be 
neglected. 

The ground state properties of Ne isotopes including the odd-A and even-even nuclei from 
the RMF+SLAP calculations with PK1 are listed in Table II, including the total binding energy 
E, the binding energy per nucleon E/A, the rms radii for neutron, proton, and matter, R n , 
R p , and R m , and the quadruple deformation for neutron, proton and matter, /3 2n , P 2p and 
/32 m - The binding energies per nucleon for Ne isotopes from the RMF+SLAP calculations are 



shown in Fig.3, in comparison with the available data 52]. It can seen that for the ground 
state, the RMF+SLAP calculations, similar as the RMF+BCS calculations [50], give very 
good description of the data, except for the proton-rich nucleus 18 Ne. However, RMF+SLAP 
approximation can do more than RMF+BCS calculations, as will be seen in the following. 
One- and two-neutron separation energies defined as, 

S n (Z, N) = E(Z, N) - E(Z, N-l), (16) 

and 

S 2n (Z, N) = E(Z, N) - E(Z, N-2), (17) 

are sensitive quantities to test a microscopic theory, where E(Z, N) is the binding energy of 
nucleus with proton number Z and neutron number N. 



In Fig. 4, the one- and two-neutron separation energies from the RMF+SLAP calculations 

n 

(filled circles) with PK1 in comparison with the data (open circles) [52] . From the figure, it 
is shown that the experimental odd-even staggering is well reproduced in the RMF+SLAP 
calculations for either the one-neutron separation energy or the two-neutron separation energy. 

In Fig. 5, the deformations /?2 for Ne isotopes from the RMF+SLAP calculations are given, 
including the data [53]. It should be mentioned that as the prolate or oblate shapes can not 
be distinguished experimentally, the assumption of shape oblate for 24 Ne has been done. First 
of all, the tendency for the shape change with neutron number in Ne isotopes is correctly 
reproduced. The RMF+SLAP calculations give even better results than the RMF+BCS calcu- 
lations [5^] . From the figure, it can be seen that the deformations are generally underestimated, 



particularly for the proton rich side. In Ref. 53], the deformation parameter (3 is determined 
from B(E2) j for the ground state to 2 + state 

P = (4tt/3ZR 2 )[B(E2) t /e 2 ] 1/2 , (18) 

where the radius -Rn has been taken to be 1.2 A 1 / 3 fm. However, it has been well demonstrated 
that the experientia! .elation of radius is on ly suitable for m edi_ and heavy nuclei Q. 
While for light nuclei the estimation for the radius Ro is too small, this leads to a large defor- 
mation (3 with Eq. (jl8|) . With this factor considered, the consistency between the calculation 
results and the deformation data will be more agreeable. 

The neutron, proton and total pairing energies from the RMF+SLAP calculations are dis- 
played in Fig. 6 as the open circles, open squares and filled circles respectively. It can be seen 
that the odd-even staggering has been clearly revealed for the neutron pairing energy due to the 
correct treatment of the Pauli blocking effects, while the proton pairing energy varies smoothly 
with the neutron number N. The large neutron pairing energy at N = 20 in present calcula- 
tions indicates disappearances of the traditional magic number at N = 20 in the Ne isotope 
chain. Although from Fig. 5 the nucleus 30 Ne is nearly spherical (fa ~ 0.1), the magic number 
iV = 20 disappears. This indicates that the disappearance of the magic number here is not due 
to the deformation, instead it might come from the fact that the neutron Fermi surface in 30 Ne 



is very close to the threshold of continuum, as noted and discussed in Refs. 



In order to justify the description for nuclei close to the drip li ne, the analytical continuation 



should be adopted 



in coupling constant (ACCC) method or the S-matrix method 
to resort the resonance states. The works along this line is in progress. 

The unified description for both the ground and excited states is another advantage of the 
RMF+SLAP. By diagonalization of the pairing Hamiltonian in the corresponding subspace 
with fixed K, parity tt, and seniority s, the spectra of multi-particle collective excitations can 



be obtained. 

As an example, the collective excitation spectra for 24 Ne are shown in Fig. 8, where the 
subfigures (a), (b), (c) and (d) respectively represent the spectra of positive and negative parity 

Q + 

for neutron, and those for proton. The calculated + excitation spectra include En 1 =5.622 

o+ 

MeV for neutron and Ep 1 =4.912 MeV for proton, which are very close to the only known 
+ 

experimental data E e l P =A.7QA MeV. From Fig. 8, it can also be seen that the RMF+SLAP can 
also be used to investigate the excited states with high-.fr value. It may provide another sel 



consistent microscopic description and prediction of the interesting high-.fr energy trap [5.7, 58]. 



IV. SUMMARY 



The Shell-model-Like APproach (SLAP) has been applied to treat the pairing correlations 
in the framework of RMF theory, which can exactly treat the Pauli blocking effects. The 
formalism and numerical techniques are presented in detail and the Ne isotopes are taken as 
an example to demonstrate the applicability of the method. 

First of all, the RMF+SLAP with effective interactions PK1 has been used to describe the 
ground state properties of the Ne isotopes, including the binding energies, one-neutron and two- 
neutron separation energies, deformation, pairing energy and occupation probability of single 
particle orbital, etc.. The calculated binding energies from the RMF+SLAP agree well with 
the data and are comparable with the RMF+BCS calculations. The odd-even staggering is well 
reproduced, as shown in the neutron pairing energy and one-neutron or two-neutron separation 
energy. The tendency for the change of the shapes with neutron number in Ne isotopes is 
correctly reproduced. Furthermore the RMF+SLAP calculations give better results than the 
RMF+BCS calculations for 26 Ne and 28 Ne. 

The other advantage of the RMF+SLAP method is the description of both the ground 
and excited states simultaneously. As examples, the collective excitation spectra for 24 Ne with 
positive and negative parity for both proton and neutron are presented. Its potential application 
for the interesting high-iT energy trap has been demonstrated. 

In conclusion, SLAP is a good method to treat the pairing correlations in the RMF theory, 
which can provide a good description for the properties in not only the ground state but also 
the excited states. 
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TABLE I: The binding energy, radii and deformation of ground state in Ne obtained from the 
RMF+SLAP* and RMF+SLAP calculations with PK1, in comparison with the RMF calculation and 
data available. Listed are total binding energy E, binding energy per nucleon E/A, neutron, proton 
and matter root mean square radii, R n , R p and R m , and the quadruple deformation parameters for 
the neutron, proton and matter distributions, fo n , fo p and fo m -, respectively. 





E E/A 


Rn Rp Rm 


fon fop fom 


EXP 

RMF 
RMF+SLAP* 
RMF+SLAP 


-191.836 -7.993 
-189.283 -7.887 
-190.723 -7.947 
-189.974 -7.916 


2.985 2.756 2.892 
2.994 2.757 2.898 
2.994 2.758 2.898 


-0.45 

-0.278 -0.238 -0.261 
-0.278 -0.236 -0.261 
-0.275 -0.235 -0.258 



TABLE II: The ground state properties of Ne isotopes obtained from the RMF+SLAP calculations 
with PK1, including the binding energy E, binding energy per nucleon E/A, neutron, proton and 
matter root mean square radii, R n , R p and R m , and the quadruple deformation parameters for the 
neutron, proton and matter distributions, fan, 02p and /?2 m , respectively, with A the mass number 
and N the neutron number. 



A N 


E 


E/A 


-*^77l 


p2n 


p2p 


P2m 


18 8 


-135.584 


-7.532 


2.555 2.951 2.782 


0.000 


0.035 


0.020 


19 9 


-144.546 


-7.608 


2.686 2.840 2.768 


0.312 


0.434 


0.377 


20 10 


-157.932 


-7.897 


2.805 2.839 2.822 


0.518 


0.533 


0.526 


21 11 


-166.982 


-7.952 


2.867 2.823 2.846 


0.513 


0.507 


0.510 


22 12 


-176.806 


-8.037 


2.928 2.813 2.876 


0.508 


0.487 


0.499 


23 13 


-182.704 


-7.944 


2.958 2.789 2.885 


0.368 


0.408 


0.386 


24 14 


-189.974 


-7.916 


2.994 2.758 2.898 


-0.275 


-0.235 


-0.258 


25 15 


-194.881 


-7.795 


3.088 2.786 2.971 


0.256 


0.326 


0.284 


26 16 


-200.902 


-7.727 


3.180 2.807 3.042 


0.277 


0.327 


0.296 


27 17 


-204.405 


-7.570 


3.253 2.821 3.100 


0.224 


0.306 


0.254 


28 18 


-209.152 


-7.470 


3.333 2.836 3.165 


0.191 


0.289 


0.226 


29 19 


-211.688 


-7.300 


3.376 2.850 3.205 


0.117 


0.234 


0.157 


30 20 


-215.708 


-7.190 


3.428 2.863 3.251 


0.055 


0.160 


0.090 
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FIG. 1: The neutron single-particle levels in Ne in RMF with PK1. The dashed line represents the 
Fermi surface A, which is also the last occupied orbit in the simple filling approximation. 
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FIG. 2: The occupation probabilities for the neutron single-particle levels in Ne and 5 Ne obtained 
from RMF+SLAP calculation with PK1, where the occupation probabilities of levels above [202]5/2 + 
for 24 Ne and those of levels above [200] 1/2+ for 25 Ne are multiplied by a factor of 10 for clarity. 
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FIG. 3: Binding energy per nucleon, E/A, for Ne isotopes as functions of neutron number N obtained 
from the RMF+SLAP calculations (filled circles) with PK1 in comparison with the data (open circles). 
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FIG. 4: One- and two- neutron separation energies, S n and S2n, for Ne isotopes as functions of neutron 
number N obtained from the RMF+SLAP calculations (filled circles) with PK1 in comparison with 
the data (open circles). 
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FIG. 5: Quadruple deformation, fe, for Ne isotopes as functions of neutron number N obtained from 
the RMF+SLAP calculations (filled circles) with PK1 in comparison with the data (open circles). 
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FIG. 6: Pairing energies, E pa i r , obtained from the RMF+SLAP calculations with PK1 for Ne isotopes 
as functions of neutron number ./V for neutron (open circles), proton (open squares) and total nuclear 
matter (filled circles). 
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FIG. 7: The excitation spectra obtained from the RMF+SLAP calculations with TMA in 24 Ne. The 
quantum number K is the eigenvalue for the third component of the total angular momentum operator 
J z . The subfigures (a), (b), (c) and (d) respectively represent the spectra of positive and negative 
parity for neutron and proton. 



